Gauge Invariance and Hall Terms in the Quasiclassical Equations of Superconductivity 
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This paper presents a careful derivation of the quasiclassical equations of superconductivity 
so that a manifest gauge invariance is retained with respect to the space-time arguments of the 
quasiclassical Green's function g. The terms responsible for the Hall effect naturally appear from the 
derivation. The equations are applicable to clean as well as dirty superconductors for an arbitrary 
external frequency much smaller than the Fermi energy. Thus, they will form a basis toward a 
complete microscopic understanding of the Ifall effect in type-Il superconductors. 
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I. INTRODUCTION 

The Hall effect in the vortex state of type-II supercon- 
ductors remains a matter of controversy after decades 
of intensive investigations. The eaxly phenomenologi- 
cal theories of Bardeen and Stepheru and Nozieres and 
VinenEl fail to account for the sign change of the Hall 
conductivity Jbund in the vortex state of a wide variety 
of materials.Era Also, a debate still continues about the 
forces acting on a single moving vortex.cl This state of af- 
fairs may be attributed partly to a lack of the established 
tractable microscopic equations with which one could test 
the validity of various phenomenological models numer- 
ically. Especially, the standard quasiclassical equations 
of superconductivity, i.e. one of the most powerful meth^ 
ods for nonequilibrium superfluids and superconductors, □ 
have been known to be unable to describe the phe- 
nomenon. 

Efforts have been made recently to include terms re- 
sponsible for the Hall effect in tic quasiclassical equa- 
tions. Larkin and Ovchinniko"vl3 incorporated higher- 
order effects arising from the particlc-holc asymmetry 
near the Fermi level, but the main terms correspond- 
ing the normal-state Hall effect are still missing in their 
equations. Kopniru obtained kinetic equations with the 
desired Hall temis and used them to discuss the flux- 
flow Hall effectJj However, their applicability is limited 
to clean superconductors with slow time variations, due 
to his transformation to Green's functions [see Eq. ( pd| ) 
below] which may not be suitable for deriving the equa- 
tions for high- frequence disturbances. More recently, 
Houghton and Vekhtert2l also reported an extension, but 
there seem to be a couple of unsatisfactory points. First, 
the obtained equations do not carry a manifest gauge 
invariance with respect to the space-time arguments of 
the quasiclassical Green's function g. The second point 
lies in their derivation process: They first define the lo- 
cal one-particle energy ^ = ^(p— |A) which depends on 
the position R as well as the momentum p through the 
spatial dependence of the vector potential A = A(R). 



Whereas they expand this ^ with respect to A to get the 
linear field dependence of the equation, the solution g to 
the equation is formally defined by the integral of the full 
Green's function G over the unexpanded ^ = ^(p— |A). 
The validity of this procedure is not entirely clear. 

With these observations, we here present a careful 
and more straightforward derivation of the quasiclassi- 
cal equations which can fully describe the Hall effect of 
the vortex states and form a firm basis for any detailed 
numerical studies. A key ingredient lies in the intro- 
duction of a new transformed Green's function [see Eq. 

below] whose gauge change can solely be expressed 
with respect to the slowly varying space-time coordinate; 
they are different from those used by Kopninu and en- 
able us to obtain completely gauge-invariant quasiclassi- 
cal equaticps. The idea goes back to the original wopk. 
of Gor'kovEl and was extensively used by Eilenbergenia 
prior to his derivation of the quasiclassical jequations.llS 
On the other hand, Larkin and OvchinnikovHil presented 
a more compact derivation of the quasiclassical equations 
with the left-right subtraction trick without recourse to 
the transformed Green's function. Those two approaches 
certainly provide the same equations at the lowest order 
of the approximation, ft turns out, however, that the ad- 
vantages of the two approaches have to be combined to 
proceed further with the gauge-invariance at every order 
of the approximation for a systematic derivation of the 
Hall terms. Indeed, the quasiclassical equations will be 
obtained here by applying the left-right subtraction trick 
to the Dyson-Gor'kov equation for the new transformed 
Green's function. 

This paper is organized as follows. Section |ll| writes 
down the Dyson-Gor'kov equation for the conventional 
retarded Nambu Green's function, followed by an intro- 
duction of a new Green's function whose gauge change 
can be expressed only with respect to the center-of-mass 
coordinate. Section III transforms the space-time deriva- 



tives of the Dyson-Gor'kov equation into an expression 
of using the new Green's function. Section IV performs 
a similar transformation to the self-energy part of the 
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Dyson-Gor'kov equation. Section ^ collects the results 
of Sees, ni and |^ to write down the Dyson-Gor'kov 
equation for the transformed Green's function, and sub- 
sequently derives the retarded quasiclassical equations by 
the left-right subtraction trick. Section ^ presents the 
exte nsio n to the advanced and the Keldysh parts. Sec- 
tion VI] concludes the paper with several remarks. 

We put h.= l throughout, and denote the light velocity, 
the electron charge, and the electron bare mass by c, 
e(< 0), and m, respectively. 



II. GREEN'S FUNCTIONS 

We will consider the elements of the Nambu-Keldysh 
matrixB separately as it turns out to be more convenient 
than handling the matrix itself. We first focus on the 
retarded part to describe the derivation, and then carry 
out the extension to the advanced and the Keldysh parts. 
We hence drop the conventional superscript ^ signifying 
"retarded" while the distinctions are unnecessary. 

Let us define a couple of retarded Green's functions by 

G„;3(l,2) = -Z0(ti-i2)(Wa(l),^]5(2)}), (1) 



F„^(l,2) = -i9ih~h){{M^),Mm . (2) 

where 1 = ri^i specifies the space-time coordinate, a/3 
are the spin indices, and {A, B} = AB+BA. To suppress 
the spin indices we introduce 2x2 matrices G and F_ by 
{G)ap = Gafi and {F_)ai3 = Faf3. Using G and F_, we next 
define a 4 x 4 Nambu matrix b}0 



G(l,2) 



G(l,2) £(1,2) 
-F*(l,2) -G*(l,2) 



and the corresponding self-energy matrix by 
E(l,2) 



S(l,2) A(l,2) 
-A*(l,2) -E*(l,2) 



(3) 



(4) 



They are changed for ~* "00(1) exp[i|x(l)] ^ 

G(l,2) ^ exp[ifx(l)T3]G(l,2)exp[-ifx(2)f3],'where 
is defined by 



T3 = 



1 

-1 



(5) 



with 1 and denoting the 2x2 unit and zero matrices, 
respectively. They satisfy the Dyson-Gor'kov equation: 



*^-e$if3)G(l,2) 



Hi 

-HI 



G(l,2) 



S(l,3)G(3,2)d3 = (5(l-2)i, 



(6) 



Here $ is the scalar potential, 1 is the 4x4 unit matrix, 
and Hi is defined by 



2m \ OTi c 



1 



(7) 



with /i the chemical potential and A the vector potential. 

We now introduce the key quantities, i.e. a couple of 
new Nambu matrices, by using a nonlocal gauge trans- 
formation as 

G(l,2) = cxp[^/(i^,rl)f3]G(l,2)expH/(i^,f2)f3], (8) 

S(l,2) =exp[i/(i?,ri)f3]S(l,2)expH^(i?,r2)f3]. (9) 

Here ri = (rfi,ri) is the four vector ,0 R = Ri2 = ^iri + 
r2), and / is defined by 



/(i?,ri) 



A{s) ■ ds, 



(10) 



where A= — A) denotes the covariant electromagnetic 
potential and ds is taken along the straight line. 

The quantities G and S defined above have a de- 
sired property that only the center-of-mass coordinate 
R is relevant in the gauge transformation ipai^) ~^ 
i/'q(1) exp[i|x(l)]. Indeed, G is changed as G(l,2) 
exp[i|x(^)f3]G(l, 2) exp[— i|x(^)'^3]- Proceeding with 
G and S, we are led to the equations with a manifest 
gauge invariance with respect, to R, as seen below. In- 
deed, Levanda and FleurovtZl has successfully derived 
normal-state kinetic equations with a manifest gauge in- 
variance by using the component G of Eq. . 

It is worth pointing out the difference of the above G 
from that used by Kopnin. As mentioned in Introduc- 
tion, KopnirJa has derived kinetic equations based on a 
couple of transformed Green's functions similar to En . 
(^. However, his phase factor is different from Eq. ( |l0[ ) 



lK{R,ri) = ~e [ ^rit')dt' + -[ A{r'ti) ■ dr\ (11) 

where dr' is along the straight line; see Eq. (18) of Ref. 

^. The present G may be more advantageous, since its 
gauge transformation property is expressible only with 
respect to R. Indeed, it will enable us a more systematic 
and comprehensive derivation of the Hall terms. 

It is convenient for later purposes to introduce the two 
functions: 

r-l 



£i{u) = I d-qe 



f2(u) = dr) dCe'^" = 



1 



1-u 



Then Eq. (|T^) can also be written as 
with f= ri2 = ri — r2 . 



(12) 



(13) 



(14) 
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III. SPACE-TIME DERIVATIVES 

Let us rewrite the first two terms on the left-hand side 
of Eq. (|) with respect to Q of Eq. (|). The following 
identities are useful for this purpose (j = 0, 1, 2, 3): 



4c 



dAj{R) dAkiR) 



dRk 



dRi 



dA,{R) dAk{R) 



dRk 



OR, 



(15) 



(16) 



Notice the differences of ±£'2(— jw) between Eqs. (|1 
and (pj|), and also between Eqs. (p^ ) and (|20|). 

We now introduce the Fourier transform of G by 



G(p£,Rr)= / G(l,2)e-'(P '"--'*)drdt 



G(pe, RT) F(pe, RT) 

-F*(-p-e,RT) -G*(-p-£,RT) 



(21) 



where the arguments e and —e both have an infinitesi- 
mal positive imaginary part. We also define the gauge- 
invariant derivatives dr and c?r by 

^ : on G, G*, S, S*, E, h 

dT={ ^ + 2ie$(i?) :on^:, A , (22) 

^-2ie$(i?) : on A* 



where u = § • and summations over the repeated 
index k (=0,1,2,3) are implied. Using the results, the 
gauge- invariant time and space derivatives of G are trans- 
formed as 



^iI{R,ri)-iI(R,r2) 



G(l,2) 

+£2(-|i^)]E(i?).rlG(l,2), 



(17) 



and 



^iI{R,ri)-iI{R,r.2) 



.9 e 
-t- A{ri, 

OTi C 



- ^ ^ l^^^'^^^^ ^ ^^^^^^ 



G(l,2) 
t a 



+•^2 ( 2 ar) 



-h(R)xr-I](R)t 

c 



G(l,2) 



(18) 



where we have neglected spatial derivatives of both the 
electric field E and the magnetic field h. Similarly, we 
have 



AI{R,ri)+iI(R,r2) 



■9 



dt 2dT 
^2(-|i^)]E(i?)-r^£(l,2) 



2aTJ '^2(297) 



(19) 



and 



JI(R,ri)+iI{R,r2} 



d i d 



-I- A ri, 

OTi C 



Z(l,2) 



ar 2 9R c 
1 



^A{R)~-^[2£,{ 



t a \ 
2 aT> 



^2 ( 2 OT ] 



t a \ 
2 ar ) 



h(i?)xr-E(i?)t 



£(1,2). 



(20) 



a 

OR 



on G, G*, S, S*, E, h 



9r= <( ^-i^A{R) :onZ, A 
^-|-i2£A(i?) : on^:*. A* 



(23) 



OR 



Then using Eqs. (p7|)-(pO|), we can write the first two 
terms on the left-hand side of Eq. (H) with respect to 
Eqs. (|2l|)-(|2^). We finally neglect terms second-order in 
c?R, E, h, and (p — Pf)^ with pp the Fermi momentum, 
since they are smaller than th£ first-order ones by an 
order of magnitude in "small" i.e. (pf'^o)^^ etc, with 
^0 the coherence length. With these procedures, the first 
two terms on the left-hand side of Eq. (^ are Fourier- 
transformed into 



-dr G 



2 

QrTsG 

2 m 



2m 



Pf 



(P - Pf) 



T3G 



--ef2(i9T4)E.0pGf3 

^^e[2fi(-iaTa,)-f2(-§5Ta,)] 

xHl.l'lhxc?, 

m 

i „ / „ „ N Pf /I 



Ea G 



-e£2{^dTd,) ^ • ( -h X + E9e ) hGh , (24) 



where dr 



5e = and £1 and £2 are defined by 
Eqs. ( p^ and (p^), respectively. Here one should keep in 
mind that dx in £1 and £2 operates only on E and h. 

With no time dependence in E and h, £1 ^ 1 and 
£2 ^ ^ so that the last four terms in Eq. (^J) reduce to 

^eE-ap(3f3G + Gf3) 


+ ^e^- fihx9p + Ea,V3G + f3Gf3). (25) 
8 TO V c / 
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These terms are absent in the conventional derivations, 
which are certainly responsible for the Hall effect of both 
the normal and the vortex states. 

Before closing the section, we compare the above result 
with the corresponding one derived by KopninE Due to 
the difference between Eqs. ( |l0|) and he obtained 

instead of Eq. (E5[) the expression: 



E -9p + — • ( ihx9p 

ni \c 



G. 



(26) 



See Eq. (23) of Ref. | In addition, whereas Eq. (H) 
are given in terms of the gauge-invariant derivatives of 
Eqs. ( p2| ) and (|2^), terms such as Vf' and dP/dt appear 
in Kopnin's Eq. (21). It should also be noted that Eq. 
( p^ ) is free from the assumption of the slow time varia- 
tions and is applicable to the cases of arbitrary external 
frequencies, as long as they are much smaller than the 
Fermi energy. 



IV. SELF-ENERGY TERMS 

We next consider the following self-energy terms ap- 
pearing on the left-hand side of Eq. ^) : 

J(l, 2) = e'^'-'^'^'^'^-''^^''-'^''^ [ E(l, 3)G(3, 2) d3 , (27) 



K{1, 2) = e*-f(^i2'^^i)-*^(«i2,r-2) / ^(1^ 3)F*(3, 2) d3 , (28) 



L(l,2)= e 



S(l,3)£(3,2)rf3, (29) 



M(l, 2) = e'^'^^'^'^^^+'^'^^'^'^^^ / A(l, 3)G*(3, 2) d3 . (30) 



Some of the main issues may be: (i) whether these terms 
can also be expressed with respect to the gauge-invariant 
derivatives of Eqs. ( p^ ) and ([2^); (ii) how the bare mass 
m in Eq. (|2^ ) is changed by the interactions; (iii) whether 
new terms arise or not besides the last four terms in Eq. 

®- ^ 

We first focus on Eq. ([27]). Writing it with respect to 

S and G, it is transformed into an expression where EG 
is multiplied by a phase factor e*'^^^^ with 



^123 : 



Ais) ■ ds 



Cv2 



dx3 



dA^ 



id: ■ (31) 



Here the contour G123 is given in Fig. 1 (a) , and we. have 
used the Stokes theorem to obtain the second line,E3 with 
the infinitesimal surface element df^'' {j,k = 0, 1, 2, 3) de- 
fined by 




Cl \ \ C 



^ (a)ZG ' ^ (b)AF* " 




(c)ZF 




(d) AG 



FIG. 1. Paths of the phase integrals 

dp' = {n2duy{Ti^dvf - {r^2du)\ri^dvy 
(0 < u < 1; 0<v <u) 

To evaluate Eq. (pl|), we expand Bjk = 
R = R12 as 

Bjk{R +{u~ i)r32 + {v- i)fi3) 
{u - i)" {v 



(32) 



dxi 



from 



oc 00 



EE 



2) 



ir32 ■ d.nns ■ d.r 



n=0 n'=0 

xB.kiR). (33) 

We then find that, with the approximation for E adopted 
below in Eq. (^7[), the terms (fia • 9^)" with n' > 1 can 
be neglected. This fact may be realized more clearly by 
noting that ri3 = (cii3,ri3) is transformed in Fourier 
space into e and p derivatives on S(p£, RT), and there 
already exists one fis in the surface element as Eq. (|3^). 
Using the approximation, the integrations over u and v 
in Eq. (^Tj) are easily performed to yield 

^123^^j[2£^{^^)^£2{'-f^)+82{~'-f^)] 



h{R) ■ (ri3 X r32) + E(i?) • (ri3i32 - ii3r32) 



(34) 

where £1 and £2 are defined by Eqs. (|l2|) and (^, re- 
spectively, and we have neglected spatial derivatives of E 
and h once again. 

We now introduce the Fourier transform of E through 

"dpde , 



S(l,3) 



(2^)4 



^^p(~2^' an." 



*22. _d_ 

2 ar 



S(p£,Ri3ri3)e^(P"- 



-etia) 



(35) 
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Then using Eq. (p4), we can transform Eq. (p7|) into 



-eti3)+i(p'-r32-e't32) 



xexp{|e[2fi(-iaTa,0 -^2(-iaT9e') +^2(1^79,0] 

X [ih(i?) . (9p X Sp. ) - E(^) • (ap5e' - 49p' )] } 
X exp [§ (9r • Op' - ^p ■ 9r' - aT9e' + 49t' )] 

xS(pe,Rr)G(pV,R'r')|R,=R,T'=T' (36) 

where 9t in 81 and £2 operates only on E and h. The lat- 
ter two exponentials are to be expanded to first order in 
<?R, E, and h, in accordance with Eq. (p^). To this en^ 
we use the following approximation for the self-energy:t2l 



E(pe, RT) « ReS(pFO, RT) + 1 



el 



— - — ) • (p-Pf) 1 + iImS(pF£, RT) , (37) 
a m J 

with a jthe renormalization factor and vp the Fermi 
velocity.tj We then neglect terms including R deriva- 
tives of ReE(pFO,RT) as well as e and R derivatives 
of ImE(pFe,Rr), since they are smaller at least by an 
order of magnitude in (ppfo)"^- We thereby obtain the 
Fourier transform of Eq. (^6|) as 

J(pe,RT) = J l{l,2)e-''^P"~''Urdt 

= E(pe, RT) o G(p£, RT) -U^^^EL). o^g 



P£) . ("ihxap + Ea.^G, (38) 



\ a m / \c 
where Ao B is defined by 

AoB = exp[^{d,dT'-dTde')]A{eT)B{e'T')\^,^^^^,^^ . 

(39) 

We next consider Eq. (^). Writing it with respect to 
A and £^ , it is expressed as A F multiplied by the phase 
factor: 

^i[I{Rl2,ri)-I{Ri2,r2)-I{Ri3,ri)-I{Ri3,r3)+I(R32,r^^ 

We add the four extra paths as the broken lines in Fig. 
1(b) and then subtract them. This factor is thereby writ- 
ten as 

Qi{<l>l+'t>2+4>3) + 2i[I{R23Ml2)-I{Ria,Rl2)] 

where (j = 1, 2, 3) is defined as Eq. ( |3l| ) with the con- 
tour Cj given in Fig. 1(b). Those 0j's may be trans- 
formed into expressions corresponding to Eq. (§3), but 



their explicit forms will not be required below; we should 
only keep in mind that they all include ^13 or ri3, as in 
the case of Eq. (|^ . 

We now introduce the Fourier transform of A in 
the same.. way as Eq. (^. We also use the following 
identity :Ej 



exp[-2i/(i?i3, ^12)] exp( 



£32 8 4_ t32 _§_ 

2 aR~^ 2 dT 



exp 



r32 . 

2 



(^-zf A(i?))+if (^ + 2ze<i>(i?))] . (40) 



With these prescriptions, Eq. ( p8| ) is transformed into 
-~ J {2^YJ (27r)4y^^3^'3" 

X e»(^i +02 +03) exp [1 (^R • dp, - dp ■ Or, - drd,, + d.dr' )] 
X A(p£, RT)£*(-p'-e', n'T')\^,^^^^,^^ , (41) 



where dx and Stl are defined by Eqs. (22) and (|2^), re- 
spectively. The latter two exponentials should be ex- 
panded to first order ift-|9R, E, and h. To this end, we 
use the approximationJlj 



A(pe,Rr) « A(pFO,Rr). 



(42) 



Compared with Eq. (p7D, the e and p — Pf expansions are 
stopped here at the lowest level; this difference originates 
from the smallness of (pf^o)~^- With Eq. (p^), all the p 
and £ derivatives on A vanish in Eq. (^) , so that we may 
put 01-1-02 + 03 = 0. Also, terms including R derivatives 
of A(pfO, RT) should be neglected as they are smaller by 
an order of magnitude in (pf^o)~^- We thereby obtain: 

^(pe,Rr) = y"£(l,2)c-*(P'-^*)drdi 
= A(pFO,RT)oT*(-p-e,RT). (43) 



Equations (|29|) and (^0|) are transformed similarly by 
using the contours of the phase integrals given in Figs. 
1(c) and 1(d), respectively. Especially for Eq. (^9|), we 
use the approximation: 



e 



A{s) ■ ds- - i A{s) ■ ds 



C2 



2£i{ 



2 dT ) 



C ( t32 d \ C I t32 d \ 

^2[—-ot) " '^2( TdT> 



-h{R) ■ (ri3 X r32) + E(^) • (ri3i32 - ii3r32) 
c 

(44) 

with C13 and G2 given in Fig. 1(c), which can be de- 
rived in the same way as Eq. (M). Notice the difference 



±f2(-^^) between Eqs. 



and 



Finally, the 



Fourier transforms of Eqs. (EOFand can be written 
in terms of the gauge- invariant derivatives of Eqs. (j2^) 
and (E3) as 
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i(pe,Rr) = y L(l,2)e-*(P'-^*)drdi 

: E(pe, RT) o F{pe, RT) - 1 (^^I _ Hi) . S^F 

2 \ a m J 



-f2(i5Ta,)]E.9pF 



V a m / 



(45) 



and 



A(pfO, RT) o G*(-p-£, RT) , 



(46) 



respectively. 

We now collect Eqs. (g|), and (g^) into a 

Nambu-matrix form with respect to G(p£,RT) and 



E(pe,RT)= / E(l,2)e-*(P '"-"*)drdt 



E(pe, RT) A(pe, RT) 

-A*(-p-e,Rr) ~E*(-p-£,Rr) 



(47) 



With the abbreviations J = J(p£, RT) and J* = 
J*(— p— £, RT), etc, the third term on the left-hand side 
of Eq. (|^) can now be written as 



-J + K -L + M 
^L*+M* -J*+K* 



-E(p£, RT) o G(p£, RT) + 1 - Hi) . OrtsG 



1 



1- - )£2{^dTd,)l]-dpGT3 



-e[2£,{-l^dTde)-£2{-^,dTde)] 



Vp 




a 


m 1 


' Vp 


_Pf 


. a 


m 



^hxap + ESe I G 



-hxSp+EaejfgGfg, 

(48) 



with 



S(P£,RT)= (l-i)£i+(^-Pl) .(p-pF)f3 

RcS(pfO, RT) A(pfO, RT) 
-A*(-pFO,Rr) -RcS(-pFO,Rr) 

ImS(pF£, RT) 

ImE(-pF-£,RT) 



(49) 



One should keep in mind that dx in £i and £2 of Eq. ( |4^ ) 
operates only on E and h. 

Several comments are in order before closing the sec- 
tion. First, Eq. (|4^ ) is expressed entirely in terms of the 
gauge-invariant derivatives of Eqs. ( p2| ) and (p^). Second, 
it exactly contains those terms which turn the bare mass 
of Eq. ( p^ ) into the effective mass m* = (dvp / dpF)~^ . 
Finally, aside from this change of the bare mass into the 
effective mass, no new terms arise in Eq. (^8|) besides the 
last four terms of Eq. (p3). 



V. QUASICLASSICAL EQUATIONS 

Adding Eqs. ( p4[ ) and ( ^8| ) yields the Fourier transform 
of the left-hand side of Eq. (|6|) . The corresponding right- 
hand side is just the unit matrix 1. Noting (£ -I- |9t)G — 
£ o G with o defined by Eq. (^) , we obtain the left-hand 
Dyson-Gor'kov equation as 



(eI - o-fs) o G - ^fsG + -VF • OrtsG 

+ '-e[2£i{-^dTd,)-£2{-^dTd,)]E ■ dpT^G 

+ '-e£2{^dTd,)E-dpGT3 



[2£i{-^dTd,) -£2{-idTd,)] 



xvF-l -hxdp + Ed, I G 
-^e£2(fM)vF-Qhx9p + E4) TaGfa 



where ^ = vf ■ (p — Pf), and a is defined by 



a(p£,RT) = af ^-Ai) i 



(50) 



ReS(pFO,Rr) A(pfO,RT) 
-A:(-pfO, RT) -ReE(-pFO, RT) 

ImS(pF£, RT) 

ImS(-pF-£,RT) 



T3 



T3 . 



(51) 



The corresponding right-hand equation may be derived 
similarly. It can also be obtained from Eq. ( pO[ ) by: 
(i) taking its Hermitian conjugate with noting the re- 
lations G^^pe.RT) = GA(pe,RT) and ±^^pe,RT) = 
S^(p£,RT), where ^ denotes "advanced"; (ii) formally 
changing ^ to ^. The result is: 



Go{ei- o-fa) - ^Gt3 - -vf • OrGts 



[2£,{^dTd,)~£2{^dTd,)]E ■ dpGh 



--ef2(-|M)E-apf3G 
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xvF-( -hx^p + E4 1 G 



: ai. 



I^T^e) VF • ( -h X ap + E I faGfg 



(52) 



Let us rewrite the above two equations with respect to 
G' = fzG. We next operate fa from the left and the right 
sides of En. (|5^) , and subtract the resulting equation 
from Eq. (|5C|). We then perform the integration over ^, 
neglecting all the ^ dependences except that of G". To 
this end, let us define the quasiclassical Green's function 

by 



g(p£,RT) = — 

flTT 



f3G(pe,Rr) cos(^0+)dC 



2(pe,Rr) /(pe,RT) 

-/*(-p-e,RT) -g*(-p-£,Rr) 



(53) 



with 0+ an infinitesimal positive constant. We also 
take the following procedures to get the final equations: 
(i) Rewrite c?p = dp^^ +vf^ with p|| the component on 
the energy surface ^. (ii) Notice Vf x c?p|| = vf x c?p and 



-(VFXh) ■ — 

c ap 



Of = 0. 



(57) 



(58) 



Without the terms with Og and O/, Eq. ( p4| ) reduces to 
the standard quasiclassical equations. By applying Eq. 
(54) to the normal state, one can show easily that Eq. 
(57) indeed describes the normal-state Hall effect. Thus, 
Eq. (^J) is expected to bring a consistent understanding 
of the Hall effect through the superconducting transition. 
Notice that the same Fermi velocity vf is relevant in both 
the acceleration term evp • and the Lorentz-force 
term |(vFxh) • ^ of Eq. (57), even after the correlation 
effects have been incorporated. 

A key point in the above derivation is that the terms 
with E • d-o,, turn out to be smaller than those with 



Vf ■ E9e by an order of magnitude in 
those terms in Eqs. (24), (pq), and ( 



(pF^o)^^- Hence 
yj) could be ne- 
glected from the beginning. Tracing back the derivation, 
it then follows that the e expansion in Eqs. (^^ and ( |4^ ) 
are unnecessary, so that a in Eq. (|54[) may have a more 



general e dependence than Eq. (|51|) as 



o-(pe, RT) = a 



2m 



l + aI](pFe,Rr)f3. (59) 



^^3G(pe, RT) cos(eO+) = . 



(iii) Neglect terms with E • cJpn , since they are smaller 
than those with vf • Ec^g by an order of magnitude in 
(PF'fo)"^- (iv) Make use of the integral expressions of 
Eqs. (|l|) and (|l|). 

Thus, the quasiclassical equations are obtained as 



(54) 



where {A, B] = AB - BA, [A, B]o = A o B - B o A, and 
Og and Of are defined by 



o„ 



d_ 
5p 



-VFxh(R,T-f779e) 
evF-E(R,r-i7?4)^|, (55) 



Of 



drj 



-VFxh(R,T- 



-FevF-E(R,T-|r/9e 



d_ 

de 



d_ 
9p 

(56) 



Hence Eq. (54) can also be used, for example, for a system 
with a strong electron-phonon interaction where there 
may be a strong e dependence in S. 

Equation ( p4|) carries manifest gauge invariance, i.e., 
it remains unchanged in the gauge transformation g — > 

exp[ifx(^)-r3]5exp[-ifx(^)T'3], and 
A-l-^. This is certainly a desired property to provide a 
support for the validity of the Pjcesent equations. Com- 
pared with the result of Kopninjj Eq. ( |5^ ) are more ad- 
vantageous in its wide applicability, i.e. it can be used for 
clean as well as dirty superconductors in arbitrary exter- 
nal frequencies much smaller than the Fermi energy. In 
addition, the terms with Og and Of are also present in 
the retarded and the advanced parts of the equations; 
those terms were neglected by Kopnin who considered 
only the static case of w = 0, but may have an important 
role in the vortex dynamics of finite external frequencies. 
One can also show that Eq. (|5j) agrees in the static limit 
to the equations obtained by Houghton and Vekhter,E3 if 
a due care is taken in the gauge choice and terms next 
order in (pfCo)~^ (i-e. terms with R and p derivatives of 
(7 and A) are neglected in their Eq. (53). Thus, Eq. ( |5^ ) 
clarifies the applicability of their Eq. (53) that it is valid 
only in the static limit; for example, the first term in the 
square bracket of Eq. (Bq) is absent in their equation. 
It follows from Eq. (M) that v = g o g satisfies 



which operate on g and / of Eq. (q3|) , respectively. If E 
and h are time independent, these operators acquire the 
simple expressions: 



[er-i - (T, z>]o + ivF • dufi 
-\{9. 05{f3,5}+0/[T-3,g]}o, 



(60) 
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with {A, B}o =AoB + BoA. In the absence of the 
right-hand terms, this equation tells us that if £> — 1 at 
some space point, as in the uniform cases, then vp • djii> 
vanishes so that i> does not change along the straight- 
line path parallel to vp; we may thereby conclude v = 1 
everywhere. However, this normalization condition no 
longer holds generally in the presence of the right-hand 
terms. This does not cause any trouble, however, and we 
only have to solve Eq. ([sj) with imposing the condition 
that z> ^ 1 as £ — > oo or E ^ 0. 



VI. EQUATIONS IN NAMBU-KELDYSH SPACE 

The above result for the retarded Green's function 
can easily be extended to the advanced and the Keldysh 
parts.El 

Let us define the advanced Green's functions by 
Gtp{l,2)^ie{t2-t,){{ij^{l),i^l{2)}) 



G';3q(2, 1) , 



FA,(l,2)^i0(t2-ii)({^„(l),^M2)}) 



-i^Ja(2,l). 



(61) 



(62) 



We also define the Fourier transform and the quasiclas- 
sical Green's function as 

G^(pe,Rr) = j d'^(l,2)e-'(P '--'*)drdt 



G'^(pe,RT) F^{pe,-RT) 
-£^*(-p-e,Rr) -G^*(-p-e,Rr) 



(63) 



and 

5^(p£,Rr) 



air 



T3G^{pe,RT)cos{^0+)d^ 



ffA(pe,Rr) /^(pe,RT) 



(64) 



respectively, where the arguments e and — e carry an in- 
finitesimal negative imaginary part. 

As for the Keldysh part, we start from the basic defi- 
nitions: 

GK^(1,2) EE -j([V'„(l),^^(2)]) = -GK*(2,1), (65) 



(1, 2) EE Mm = 'Fpai^, 1) . 

We then introduce the Nambu matrix by 

G^(l,2) Z^(l,2)' 
£'^*(1,2) G^*(l,2) 

Its Fourier transform is defined by 



G^(l,2)^ 



(66) 



(67) 



G^(pe,Rr)= G (l,2)e-*(P'-^*)drdt 



G^(p£,RT) F'^(p£,RT) 
£'^*(-p-£,RT) G^*(-p-£,RT) 

and the quasiclassical Green's function by 



(68) 



5^(p£,RT) 



an 



f3G^(p£,Rr) cos(^0+)dC 



5^(p£,Rr) /^(P£,RT) 
/^*(-p-£,Rr) g^*{-p-e,RT) 



(69) 



The Keldysh self-energy matrices and ct^ are defined 
similarly as Eqs. ( |68| ) and (|6|), respectively. 

We now introduce as usual three Keldysh matrices by 



(70) 



Then the equations for , g^, and g^ can be put into a 
compact form as 













6 ' 


6 g^_ 


, a = 


6 


, T3 = 


6 





+ |c?9{^3,.9} + ^0/[r3,5] =6, 



(71) 



where Og and Of are defined as Eqs. ( pq ) and (jSq), re- 
spectively. 

Those quasiclassical Green's functions satisfy 

[g^(p£, RT)]t = -f3.9^(p£, Rr)f3 , (72) 

[5^(P£, RT)]t = f35^(p£, RT)f3 , (73) 

[gi^(p£, RT)]T = f2gA(-p-£, Rr)f2 , (74) 

[5^(p£,RT)]T = f2.9K(-p-£,Rr)f2, (75) 

with denoting the transpose. These relations originate 
from [GR(p£,Rr)]t = G^[pe,RT), [G^{pe,'KT)\^ = 
-GK(p£,Rr), [Gi^(p£,RT)]T ^ -fiGA(-p-£,Rr)fi, 
and [GK(p£,RT)]T = -fiG'^(-p - £, RT)fi, respec- 
tively. Similarly, the self-energy matrices are shown to 
have the symmetry: 

[a^(p£, Rr)]t = f3(T^(p£, Rr)f3 , (76) 

[aK(p£, Rr)]t = -f3aK(p£, RT)f3 , (77) 

[a'^(p£, RT)]T = T2a^{-p-e, RT)f2 , (78) 

[aK(p£, RT)]T ^ f2aK(-p-£, Rr)f2 . (79) 
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VII. CONCLUDING REMARKS 

We have presented a systematic derivation of the qua- 
siclassical equations based a nonlocal-gauge-transformed 
Green's function (||). This enabled us to retain the gauge 
invariance in terms of the center-of-mass coordinate R = 
(RT) at every stage throughout the derivation. Equation 
(H) with Eqs. (H), fg]), (|9|), (H), and is the 

main result of this paper. It naturally carries a manifest 
gauge invariance, i.e., it remains unchanged in the gauge 
transformation 5— »exp[i|x(^)r3]^exp[— i|x(^)''^3]; * 
andA 



A + ^ . This is certainly a desired 



property to provide a strong support for the validity of 
the present equations. Also, the terms responsible for 
the Hall effect are automatically present in the operators 



Og and Of. Indeed, by applying Eq. (54) to the nor- 



mal state, one recovers the normal-state Hall effect. It 
should also be noted that Eq. ( p^ is applicable to band 
electrons; this may be shown by using in the derivation 
the anisotropic self-energy S(pe, RT) where the effect of 
the periodic lattice potential is incorporated. 

Compared-^ith the results of KopninQ and Houghton 
and Vekhteill3 which are valid only in the static limit, 
as discussed in the paragraph below Eq. (|59|), Eq. (|4|) is 
more advantageous in its wide applicability that it can be 
used for clean as well as dirty superconductors up to the 
external frequencies comparable with the energy gap. In 
addition, terms with Og and Of are also present in the 
retarded and the advanced parts of the equations; those 
terms were neglected by Kopnin who considered only the 
static limit of w = 0, but may have an important role in 
the cases of finite external frequencies. 

Thus, we have derived an equation which forms a firm 
basis for detailed studies of the Hall effect in the vortex 
states. Solving Eq. ( |5^ ) is expected to bring a compre- 
hensive understanding of the Hall effect in type-II super- 
conductors. 
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